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Q ■ Abstract. We study the phenomenology of hybrid scenarios of neutrino dark energy, 

where in addition to a so-called Mass Varying Neutrino (MaVaN) sector a cosmological 
constant (from a false vacuum) is driving the accelerated expansion of the universe 
today. For general power law potentials we calculate the effective equation of state 
parameter Weff{z) in terms of the neutrino mass scale. Due to the interaction of the 
^ ' dark energy field ( "acceleron" ) with the neutrino sector, Weff (z) is predicted to become 

^D ! smaller than —1 for z > 0, which could be tested in future cosmological observations. 

^L} ' For the considered scenarios, the neutrino mass scale additionally determines which 

^^ . fraction of the dark energy is dynamical, and which originates from the "cosmological 

constant like" vacuum energy of the false vacuum. On the other hand, the field value of 
("^ • the "acceleron" field today as well as the masses of the right-handed neutrinos, which 

QO . appear in the seesaw-type mechanism for small neutrino masses, are not fixed. This, in 

principle, allows to realise hybrid scenarios of neutrino dark energy with a "high-scale" 
seesaw where the right-handed neutrino masses are close to the GUT scale. We also 
comment on how MaVaN Hybrid Scenarios with "high-scale" seesaw might help to 
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$H ' resolve stability problems of dark energy models with non-rclativistic neutrinos. 



1. Introduction 

The evidence for tlie existence of a dark sector in the universe has made the present 
era of cosmology fascinating and challenging. At present, the microscopic natures of 
dark energy and dark matter are still an open question, with both components only 
probed gravitationally [H [21 [3l H]. Regarding dark matter, there are various particle 
physics candidates which may belong to the class of Weakly Interacting Massive Particles 
(WIMPS) or which may interact only gravitationally. On the other hand, for particle 
physics explanations of the observed dark energy, the main fundamental question is 
whether dark energy is a cosmological constant or a dynamical field. While in the 
former case the equation of state parameter w of dark energy is constant and equal 
to —1, in the latter case it is a function of redshift z and in general differs from —1, 
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providing a way to distinguish both fundamental types of dark energy experimentally. 
Such a deviation of w from —1 will be searched for in various future surveys, studying 
supernovae, baryonic acoustic oscillations, weak gravitational lensing effects, galaxy 
clusters, or other techniques [5]. One obstacle to relate dark energy more explicitly to 
particle physics is the fact that the amount of the observed energy density p^^ ~ 0.003 
eV is much smaller than the "theoretical expectation" . In this paper we will not address 
this "cosmological constant" problem, i.e. the suppression/cancellation of the various 
generic contributions to dark energy which are too large, but assume that it is resolved 
by some other mechanism. 

Our work is motivated by the intriguing observation that the energy density of 
dark energy is close to another very small scale in particle physics, namely the one of 
neutrino masses. The discovery of flavour conversion of neutrinos from various sources, 
interpreted within the framework of neutrino oscillations, points to two mass eigenvalues 
of the light neutrinos above about 0.01 eV and 0.05 eV [6|, while searches for neutrino 
masses from Tritium /3-decay and neutrinoless double /?-decay yield an upper bound for 
each mass eigenvalue of roughly 0.5 eV [6J. The proximity of the scale of dark energy to 
that of neutrino masses has inspired the proposal of the so-called Mass Varying Neutrino 
(MaVaN) scenario [7|. In this scenario, the relic neutrino density contributes to the 
effective potential for the dynamical field which is called "acceleron" . One consequence 
is that neutrino masses vary on cosmological time scales and, since the relation between 
neutrino masses and pj^^ remains valid in earlier cosmic epochs, also dark energy behaves 
differently from a cosmological constant. Compared to other dynamical models of dark 
energy (e.g. quintessence [8l[9l[T0]), in the MaVaN scenario one does not have to choose 
a mass for the dynamical field of the order of the present Hubble parameter Hq. In fact it 
can be much larger, even of the order of the neutrino mass scale, due to the stabilising 
effect of the contribution of the relic neutrino density to the potential. As a result, 
the field adiabatically tracks the minimum of the effective potential whose evolution is 
controlled by the time evolution of neutrino number density. 

The MaVaN models have been investigated in many studies, e.g. regarding possible 
experimental signatures and constraints [HI [T2|, [131 [HI [13 UHl [13 [HI [HI [20] as well as 
model-building issues [211 1221 ESI [25l [26]. It has turned out that while there are many 
interesting possible signatures and attractive features, the scenario is tightly constrained 
by the requirement of consistency with late time structure formation. In particular, it 
has been pointed out [27] that if neutrinos are non-relativistic today, due to fifth force 
effects the neutrinos would cluster at late time and finally form "neutrino nuggets" 
which would spoil the dark energy behavior of the neutrino fluid. 

The goal of this paper is to investigate the phenomenology of a generalised scenario 
which we will refer to as MaVaN Hybrid Scenario, where in addition to a MaVaN sector 
a cosmological constant (from a false vacuum) is driving the accelerated expansion of 
the universe today. As we will show, for a generalised power law potential in the MaVaN 
sector the effective equation of state parameter Weff{z) as well as the fraction to which 
the dark energy is of dynamical nature is determined by the neutrino mass scale. We will 
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also analyse the possibility to realise neutrino dark energy with a "high-scale" seesaw 
mechanism, where the right-handed neutrino masses are close to the GUT scale, and 
comment on how such models might help to suppress the formation of "neutrino nuggets" 
and resolve stability problems of MaVaN models with non-relativistic neutrinos. For the 
main part of the paper we will focus on neutrinos that are non-relativistic today. 

The paper is organised as follows: In Section 2 we will review the basics of the 
MaVaN scenario. Section 3 contains the definition of the phenomenological framework 
and the analysis is performed in Section 4. In Section 5 we discuss the possibility of 
realising neutrino dark energy with a "high-scale" seesaw mechanism and in Section 6 we 
comment on how this might help to solve stability issues for non-relativistic neutrinos. 
Section 7 contains our Summary and Conclusions. 

2. The Original MaVaN Scenario 

In the basic MaVaN scenario a singlet fermion N (right-handed neutrino) and a 
dynamical real scalar field A (called "acceleron" ) are introduced. Their dynamics gives 
rise to the acceleration of the universe and to neutrino masses which vary with time. 
The basic form of the Lagrangian is given by 

£ D mouN + -kANN + H.c. + V{A) , (1) 

where rriD is neutrino Dirac mass, k is a dimensionless coupling in dark sector and V{A) 
is the potential energy of the acceleron field. The Standard Model neutrino u as well as 
the singlet neutrino N are both written as two-component left-chiral Weyl spinor fields. 
If kA ^ rri/), i.e. if the Majorana mass of the right-handed neutrino is much larger than 
the neutrino Dirac mass, one can effectively integrate out A^ from the low energy theory 
which leads to the effective Lagrangian 

1 TD 

CD--^uu + }i.c. + V{A). (2) 

2 

The neutrino mass m^, = |^| depends on the dynamics of A. Taking into account the 
masses of the non-relativistic relic neutrinos we now include a term proportional to the 
relic neutrino number density n^ to the effective potential. 



Veff = n. 



m 



D 



The effective potential Veff for the acceleron A has a minimum because the first 
term pushes A to larger values whereas as the second term V{A) is assumed to favour 
smaller values of A (e.g. if V{A) oc log(74/A), or oc A"). Typically, in this type of models 
a sub-eV scale is introduced in order to realize the desired values of V{A) (and finally 
of K//(^))lil As a consequence of the two contributions to the potential and since n^ 
decreases with time, the field value of the acceleron field A becomes smaller with time 

I A scenario of quintessence cosmology with growing matter component (which may be apphed, for 
instance, to neutrinos) has recently been proposed without explicitly invoking any sub-eV scale [2H1 ES] • 
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when it adiabatically tracks the ininiina of the potential (if the mass of A is above the 
Hubble parameter). The combined A and neutrino fluid drives the accelerated expansion 
of the universe. 

We remark that unlike in conventional quintessence models the scalar field does 
not roll down an extremely fiat potential but it rather gets stabilised by the finite 
neutrino density. For this reason, the slope of the potential can be much higher than 
in quintessence scenarios and the mass of the scalar field A can even be as high as the 
neutrino mass scale, whereas in quintessence models it has to be as tiny as the Hubble 
scale. 

Although the model can successfully explain the present acceleration of the universe 
with masses of the dynamical scalar field much larger than the Hubble scale (e.g. in the 
sub-eV range) and provide a connection between dark energy and the physics of neutrino 
masses, it has also some shortcomings: For example, the potential is not stable under 
radiative corrections, unless new states appear with sub-eV masses, and it still has 
to be rather fiat in order to be consistent with present constraints on the equation of 
state parameter. Furthermore, the MaVaN scenario with non-relativistic neutrinos has 
a problem with respect to the consistency with late time structure formation, as has 
been pointed out in Ref. [27]. If neutrinos are non-relativistic today, due to fifth force 
effects the neutrinos could cluster at late time and finally form "neutrino nuggets" which 
would spoil the dark energy behavior of the neutrino fluid. 

However, it has been demonstrated that under certain conditions these problems 
might be overcome [2Tj. For example, in supersymmetric theories of neutrino dark 
energy the potential can be stable under radiative corrections (for example in gauge 
mediated supersymmetry breaking where supersymmetry breaking is only transmitted 
to the dark sector via gravity such that the dark sector superpartners are light) and 
stability problems with respect to neutrino density perturbations can be avoided if 
neutrinos are highly relativistic, although then the connection to neutrino masses is 
weakened since the two observed mass scales correspond to non-relativistic neutrinos. 

3. Generalised Framework: The MaVaN Hybrid Scenario 

After having reviewed the basics of the original MaVaN scenario in the last section, 
we will now set up the framework for our analysis. It consists of a generalised MaVaN 
potential of power law type plus a constant vacuum energy contribution. In the following 
we will specify it explicitly and comment on how it may be realised in classes of models 
as well as on the connection to the (in principle) experimentally measurable neutrino 
mass scale. 

3.1. Generalised Potential for Neutrino Dark Energy 

In the following we will consider a rather general setup where neutrino dark energy 
is realised in a type of "hybrid dark energy" setup with the theory being in a false 
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vacuum at present time. This implies that in addition to the time varying contribution 
of the acceleron field there is an additional contribution of the energy density of the 
false vacuumlj Motivated by the notion of providing a rather general framework for our 
study, we will therefore consider an effective potential for neutrino dark energy of the 
form 

ViAUf = V{A) + P^i^^j +Vo, (4) 

where 

V{A) = M^-"A". (5) 

A is the acceleron fielctjjj, pi is the energy density of the SM neutrinos today and a is 
the scale factor. The main ingredients are a constant term Vq (from a false vacuum), a 
potential for A which is monotonically increasing with \A\ (in our case a general power 
law potential oc A") and an effective term UyTny, where we assume my oc 1/A^ in order 
to generalise the conventional seesaw relation. We note that the potential of Eq. ([5]) is 
not the most general form of a potential that may arise in a realistic model. Instead of 
the monomial potential that we are considering, it can for instance also be a polynomial 
or logarithmic function of the field A, or a combination of both. 

3.2. Possible Origin of the Generalised Potential 

To motivate the appearance of the constant vacuum energy in the "hybrid scenario" 
by an example, we may consider the following superpotential (after EW symmetry 
breaking) , 

W = \A{N^ - vl) + yHlvN + Kn\A^ , (6) 

where hats indicate superfields. A^ denotes the SM singlet (right-handed neutrino) 
superfield(s), and the most relevant contributions to the scalar potential stem from the 
F-terms |-FaP + l-^AfP = I-^'^atP + |2Ay4A^ + m^A + ...p. The first term contributes 
the constant vacuum energy while the second term gives a mass term (~ |Ap) to the 
acceleron field. Together with the additional effective potential term Uyrrii, ~ nj,f^/|y4|, 
where Vu is the vev of the Higgs field which couples to the neutrinos, the potential for the 
acceleron has a minimum for a non-zero field value of A. |F/vp also gives a time- varying 
mass to the right-handed neutrino, rriN ~ |^|. 

While the above superpotential may serve as a simple example which could realise 
the generalised potential of Eq. (jl]) with a = 2 and [3 = 1, it should be kept in mind 
that that there is in principle a large variety of possible models which can give rise to 
other values of a and j3 (as well as to other non-power-law forms of the potential). For 
example, different powers of A may appear in the superpotential or the neutrino masses 
may originate from a different type of seesaw mechanism (e.g. from a so-called double 

§ One such example model is the SUSY version of the MaVaN scenario proposed in [21] . 
II For simplicity we will take A as a real scalar field as in Section [2l unless stated otherwise. 
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[20] or multiple seesaw) such that n^rrii, ~ n^v'^/\A\'^, and so on. Furthermore, one may 
anticipate that the parameters a and (3 in realistic values should be positive integers. 
However, it is also possible that fractional values of the exponents appear, for example 
if the kinetic terms in a model are non-canonical and therefore a field transformation 
has to be performed in order to bring the kinetic terms back to canonical form. We will 
discuss the possible origin of the generalised potential in more detail in the Appendix. 

3.3. Connection to Neutrino Masses 

From neutrino oscillation experiments, two so-called mass squared differences are known: 
Ar/ig^ = mg — ml is known mainly from atmospheric neutrinos and Am^i = fn^ — m\ 
from solar neutrino data and from the KamLAND experiment. The experimental values 
are about Amgi = 7.6 x 10"^ eV^, jAml^^l = 2.5 x lO^^eV^ [6]. With this nomenclature 
we also know that m2 > mi, but m^ can be the heaviest or the lightest of the mass 
eigenvalues. For the sum m,^ defined as 

3 

rh^ = '^mi (7) 

we know that there is a lower bound from oscillations (saturated iimi ^ and m2 ~ 0.01 
eV, m^ ^0.05 eV) and also an upper bound from neutrinoless double f3 experiments of 
about 1.5 eV (with mi ^ m2 ~ m^ ^ 0.5 eV) [6]. Finally, we have approximately 

m^ G [0.06 eV, 1.5 eV] . (8) 

On the other hand, the lightest of the neutrino mass eigenvalues {nii or ms) can be 
arbitrarily small. 

One important question is whether this sum m,y is relevant for the MaVaN setup, 
or whether the mass of only one of the light neutrinos matters. To answer this question, 
we have to go to the full three family scenario where we distinguish two possible cases. 

In one possible situation (which we may call case I) m,^ is indeed relevant. This 
case applies when the vev of the field A gives masses to all the right-handed neutrinos 
Ni such that the part of the Lagrangian relevant for neutrino masses reads for instance 
{Xij is a coupling matrix) 

C D ^X,,AN,N, + {Y,)^,u^HlNi . (9) 

The seesaw suppressed mass matrix of the light neutrinos is then given by 

{^v)<.ti = J^ ■ (10) 

If we diagonalise m^, the effective potential reads 

Veff = J^n^^mf s),, = n,5^(mf s),„ = n,m, , (11) 

a a 

where we have assumed that n^^ = n^ for all generations a of neutrinos. 
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On the other hand, we could have the situation (which we may call case II) that 
each Ni gets its mass from a different scalar field, e.g. from the vevs of fields A, B and 
C and that we have a Lagrangian of the form 

+ Y,(y,>^HlN, ' (12) 

i 

where Xi are coupling constants. Now the light neutrino mass matrix would be given by 

^"^^)"^ = -^;a- + ^^ + ^^ • J'') 

In this example, only mi would be relevant for the effective potential for A\^ If 
consistency of the MaVaN scenario (for a particular form of the potential) would require 
a very light neutrino mass scale ^ 0.01 eV, this could be realised in case II. The relevant 
mass would then be mi for a so-called normal hierarchical spectrum or m^, for an inverted 
hierarchical spectrum. 

In the following, we will mainly be interested in case I, since it has a particularly 
close link to the observed neutrino masses which manifests itself in the upper and lower 
experimental bounds on today's value of rhu given in Eq. ([8]). We will therefore present 
our results using the notation for case I where the relevant neutrino mass scale is rh^. 
Most of the formulae can readily be adopted to case II, as long as neutrinos remain 
non-relativistic. 

4. Phenomenology of the MaVaN Hybrid Scenario 

In this section we will investigate the dynamics and phenomenology of the MaVaN 
Hybrid Scenario with a general power law potential 

V{A) = M^-^A"" . (14) 

The observed small neutrino masses are generated by a version of the seesaw mechanism, 
and we will assume a general dependence of the neutrino masses on A of the form 

mM) ~ -jp , (15) 

where v is the vev of the SM Higgs field (which is proportional to the Dirac masses of 
the neutrinos). The notation fhy is introduced in section [3731 

4.1. Background Evolution 

As a first step of our analysis we consider the background dynamics of the scalar field A 
in the MaVaN Hybrid Scenario defined in section [3l which leads to an effective potential 
for A of the form 

V{AUj = M^-A- + ^^ + Vo , (16) 



% Of course, analogously, one can realise the case where any of the mass eigenvalues of the light 
neutrinos is relevant for the effective potential for A. 
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where nl rhl = pi and where the index "(0)" indicates todays values of the 
parameters. The second term in the potential originates from the coupling between 
neutrino and the scalar field in the effective Lagrangian. In the following discussions, 
if we do not mention the values of a and /3, we will assume a = 2 and /3 = 1 
which corresponds to the "standard" MaVaN scenario, i.e. where the potential for A 
is quadratic and where the seesaw mechanism is of type I. We will also focus on the case 
where the neutrinos are non-relativistic today. 

The Friedmann equation in our scenario takes the form 

r2A/r2 _ rCDM , t^Baryons , t/Z/IN , /^^ / ^ 1 i A'2 



n,H^Ml, = '-^^ + '^^^^^ + V{A) + '-W-^\ +-A' + Vo (17) 

a-^ a-^ a-^ \ A J 2 

and the scalar field A obeys the following modified Klein Gordon equation 

We note that while the dark matter density pcdm redshifts as a~^, the time dependence 
of the neutrino energy density is non-trivial due to the time- dependence of A. 

Now we use adiabaticity, i.e. the feature of MaVaN scenario that the acceleron field 
adiabatically follows the minimum of the effective potential since long time back in the 
cosmic history. Consequently, the condition —^^ = determines the evolution of the 
scalar field to be 

A = (^P^M-^'-A "^' . (19) 

The present value of Aq can be found from the above equation and is given by 

A, = (f^M-^'-A , (20) 

which allows to derive the relation 

^ = il + z)^ , (21) 

where a = 1/(1 + z) has been used with z being the redshift. From this expression we 
can calculate the kinetic energy of the field compared to the total energy density carried 
by the field, 

iVH.,5(JL)V^)^:..,^. p2, 

pA \MpiJ \a + (3J 

We can see that if the field value of A today is smaller than the Planck scale Mpi, 
the kinetic energy of the field can consistently be neglected for the relevant time of 
the evolution of the dark energjtj. In addition, a constant vacuum energy contribution 
further reduces the significance of the kinetic energy compared to the Hubble constant 

+ On the other hand we wiU see that a high value of ^o ^ Mpi is typically required with respect to 
the stability issues for non-relativistic neutrinos (c.f. discussion in section [6]). 
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during the evolution. We will discuss the adiabaticity condition of the acceleron field in 
Section [51 

Using Eq. ( l2Ti) and comparing the potential energy V{A) to the acceleron-dependent 
neutrino fiuid density in Eq. (fT6|) . we find that the ratio of the cosmological density 
parameters of these two components is given by 

where we have substituted 3H^n^ = ^{^^ and SH^^a = V{A). We would like 
to highlight the result that the ratio of the energy densities from the interaction term 
and from the potential term V{A) is constant in time. This feature, which is basically 
a by-product of the adiabaticity condition, will be useful to derive constraints on the 
model parameters from cosmological observations and it will lead to a particularly close 
connection to the neutrino mass scale in the MaVaN Hybrid scenario with power law 
potentials. 

4.2. Dark Energy Effective Equation of State in the General MaVaN Hybrid Scenario 

Using the above results for the background evolution, we will now investigate how the 
parameters a, (3 and Vq of the model are constrained by cosmological observations, in 
particularly regarding the dark energy equation of state, as well as by the present bounds 
on (and eventually by a future measurement of) the neutrino mass scale mj,o 

We will start by deriving the effective dark energy equation of state parameter 
Weff{z). The main point here is that when the dark energy equation of state is extracted 
from cosmological data sensitive to the Hubble scale H{z), it is always assumed that 
the dark matter component redshifts as 1/a^. Massive neutrinos also contribute a small 
fraction to the dark matter today. The total dark matter density today is then given by 

(0) _ (0) , U(0) (r.A\ 

PdM — PCDM + PdM ; l^^J 

where Pcj:,^ is the amount of cold dark matter and where p^jv/ ^^ calculated from the 
actual value of the neutrino mass ml (which may be determined more precisely by 
future experiments), i.e. p^^ = pi'. However, as we have already mentioned, the 
neutrino masses vary with time and therefore this contribution does not redshift as 
1/a^. Consequently, an observer will not extract the intrinsic dark energy equation of 
state, but rather an effective one, We/f- 

More explicitly, the effective dark energy density (which defines the effective 
equation of state Wefj) is extracted from the Hubble equation (assuming for simplicity 

* We would like to note that only the terrestrial bound on neutrino masses can be applied directly 
to MaVaN models. Cosmological bounds on neutrino masses are generically relaxed for mass varying 
neutrinos (see e.g. [32j ) since the neutrino masses decrease with increasing z, as we will discuss below. 
In fact, one "smoking gun" signal of MaVaN scenarios would be a "cosmological upper bound" on the 
sum of neutrino masses (if they are assumed as constant in time) which is incompatible with todays 
lower bound from terrestrial experiments. 
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a flat universe) 

(0) (0) 1/(0) 

SH'M'pi = ^^ + ^^^^^ + ^ + p^^i , (25) 

a-^ a-^ a'^ 

where it is (wrongly, but for a model-independent analysis unavoidably) assumed that 

the neutrino density redshifts like ordinary matter. 

We would like to add a few remarks to clarify the notation: p^^| (i.e. the value of 
the quantity p^^^ today) is same as pi. . However, in an earlier epoch, the true neutrino 
energy density pi, is not equal to the quantity p^j^,^ = p^j^^j/a^. 

To calculate p^^, following [21], we rewrite Eq. flT7|) as 

(0) (0) {0)u (0)1/ 

0-2 11^2 _ PCDM _|_ f^Baryons PpM _ PpM 

a^ a^ a^ a^ 

+ '^^ + nA) + lA^ + Vo. (26) 

Comparing this relation with Eq. (!25|) yields 

"- - (^^ - #) + yW . \a^ + V^„ . (27) 

From p^^ the observed equation of state We// can be calculated via the continuity 
equation 

^ = -3iJ(l + ^e,,)p^^i. (28) 

Taking the time derivative of Eq. 0271) and using Eq. (ITS]) and the expression for p'^, we 
find 

^ = -3^^ f ^ - 1 V 3i/i^ . (29) 

dt a^ V mi°) ; ^ ^ 

We compare the above equation with Eq. ( l28l) to obtain 

^^^^ = — (0) / - ,A. \ 77 • (30) 

This is the general result for the effective equation of state parameter Weff{z) for general 
A-dependent neutrino masses m^{A) and for a general acceleron potential V{A). 

We now discuss some of the characteristic feature of Weff'- To start with, it can 
be seen from the above expression that at present time {z = 0), the pL -dependent 
term disappears from the denominator. Furthermore, in the absence of the constant 
vacuum energy term Vq, the equation of state parameter at 2; = is simply given by the 
conventional equation of state parameter for the scalar field A, 

AV2 - V{A) 

wa = -■ • -jl) 

A''/2 + V{A) ^ ' 

Since the slow-roll approximation can be applied for field values less than Mpi and since 
the kinetic energy ~ A^ is much smaller than V{A) in this case, wa is very close to -1. 
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Regarding the effective equation of state parameter Wefj for z = 0, adding tlie constant 
vacuum energy term Vq even makes its value closer to —1. 

For 2; > 0, however, the neutrino density term in the denominator of Eq. ( l30l) is non- 
zero. Furthermore, noting that in the MaVaN-type models of interest Th,^{A) < rhi , 
the p;y-dependent term is negative which implies that Weff becomes less than —1 for 
z>0. 

In summary, the effective equation of state parameter Weff{z) is defined by Eq. fl28|) . 
with p^g related to the Hubble parameter H{z) by Eq. ( l25l) . By definition, Weff{z) 
contains information about whether the effective dark energy density is dynamical or 
behaves like a cosmological constant at a given redshift z. In order to extract p^^ from 
H{z), additional information regarding the matter components has to be used. For 
example, SN data provides some information on H{z). However, Pcj:,^ + Psaryons i^ 
required as an input, in addition to the neutrino energy density p}^]^ which depends on 
the (in principle) measurable neutrino mass scale today. One possibility to obtain P(j\)}^ 
and PBaryons ^^ ^o extract it from the CMB, where the neutrinos were relativistic and do 
not contribute to the matter density. On the other hand, if information on the matter 
density from other observations, e.g. from weak lensing, is used it may already contain 
p'" if the neutrinos are clustered on the relevant scales. We note that in order to confront 
the MaVaN Hybrid Scenarios with the whole variety of (future) data from cosmological 
observations, such issues have to be taken into account carefully. Finally, in order to 
test a particular dark energy model of this type using all available cosmological data, 
one would have to go beyond the analysis of Weff{z) and fit all model parameters to the 
data. 

4-3. Dark Energy Effective Equation of State in the MaVaN Hybrid Scenario with 
Power Law Potentials 

We will now quantify the above general statements about Weff{z) for the more explicit 
scenario where the MaVaN potential has power law form, and investigate the further 
phenomenological consequences. We start with rewriting the general form of the effective 
equation of state in terms of the parameters of the MaVaN Hybrid Scenario, which (at 
this stage) are a,/3, and Vq (assuming further a fixed/measured value of rh,y). To do so, 
we use the results from our adiabatic background solution p,y = ^pA (see Eq. [23]) and 
the relation 

p^(2) =pW(l + 2)^ or m^(z) = m[,°)(l + z)"^ , (32) 

which is plotted in Fig. ([T]) for illustration. 

As the neutrino mass scale increases with time, neutrinos were relativistic at some 
earlier epoch. For example, for ^ = 2, this happens around redshift z = 100. We note 
that at this (or larger) redshift, the above-derived formulae do no longer apply since 
we have assumed non-relativistic neutrinos throughout. The field A remains stabilised 
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when the neutrinos are relativistiqj, however at values which in general differ from those 
given in Eq. flT9|) . During CMB formation time where z ~ 1000, neutrinos were highly 
relativistic and furthermore their mass was much smaller than today, which implies 
that the effects of their mass density can be neglected during this epoch. Details about 
cosmological bounds from CMB data and large scale structure observations (including 
neutrino masses) can be found, e.g., in [T9l [32] . 

Using further the constraint from the observed dark energy density today, 

(33) 



(34) 



(35) 

for the fractional contribution of Vq to the total observed dark energy density today. 
After some straightforward calculations, we finally get the expression for the equation 
of state in terms of the parameters a, /3, 7 and as a function of redshift z, 

I ^ 7 + (l-7)(l + z)^ 

Weff[z) = 3^; . (36) 

7 + (1 + f )(1 - 7)(1 + z)—^ - f (1 - 7)(1 + zf 

We can see from the above expression that for 7 = 1, Weff{z) = — 1 as this should be the 
case for a constant energy density. However for 7 7^ 1 we find that Weff{z) deviates from 
—1 for z > 0. The amount of deviation depends on the choice of parameters a, (3, 7. 
We can also see that the value of the mi does not appear in Eq. fl36|) . As we will 







Pde 


= K) + 


Pa 5 


we find that 












pZ 


-> 


-7), 


where 


we have defined the parameter 






7 = 


Pde 





show in the next subsection, fni is determined in this setup by the parameters a,/3 
and 7. In fact, we use this to trade the parameter 7 for the measurable quantity ffii, 
and re-express Weff{z) as a function of mi ' . 

4.4- The Effective Dark Energy Equation of State and the Neutrino Mass Scale 

As mentioned above, the neutrino mass scale mi can be calculated as a function of the 
model parameters a,f3 and 7. For non-relativistic neutrinos (corresponding roughly to 
TRi > 10~^ eV for i = 1,2, 3) the neutrino energy density is pi = ffii ni . As we have 
mentioned in section 13. 3[ we are considering case I where the acceleron field provides 
the masses of all three right-handed neutrinos. Furthermore, the observed dark energy 
density today is given by p\)^ = Vq + V{Aq). Using these two equations one can show 
that mi can be expressed in terms of 7 as 

(0) 

-i'°^ = (l-7)|^- (37) 

P nl 

tt See e.g. [21] for a discussfon in the original MaVaN scenario which can readily be generalised to the 
hybrid case. 
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Figure 1. The left panel shows the neutrino mass scale ffi^ — X]i=i"^«' divided by 
their present value, as a function of redshift and of a/p. The right panel shows the 
value of 7 = Vq/ pj^^ as a function of the neutrino mass scale today and of a/p. 



Plugging in the observed value p]-,^ ~ 4 x 10 ^^ eV^ and the (standard) theoretical 
prediction for the neutrino number density today, ni ~ 8.8 x 10~^^ eV^ [22l [33], we 
find 



m?) 



45(1 



oO|eV, 



(3J 



This simple expression relates the present value of the neutrino mass to the quantity 7 
which specifies which fraction of dark energy is constant vacuum energy. For fixed 7, 
the neutrino mass scale tjiI is predicted in this scenario. 

One immediate consequence of the above relation is that for 7 = 0, i.e. without 
a constant energy density contribution Vq, the predicted neutrino mass violates the 
the present experimental bounds, unless either a is very small or (3 is very large. In 
particular, for the standard case a = 2 and /? = 1, and for rfii of the order 1 eV, Vq 
must contribute ~ 99% of todays dark energy of the total observed dark energy density 
(c.f. Fig. [1]). On the other hand, for smaller a/ (3 the value of 7 decreases and a larger 
fraction of dark energy is dynamical. For example, for P = l,a = 0.1 and rhl ~ 1 eV, 
the dynamical component amounts about 22% of the total dark energy density. 

Using Eq. (138!) we can express the effective equation of state parameter of Eq. (!36ll 
in terms of ml by eliminating the parameter 7, which yields 



Weff 



1-^ 



. (0) 



3a 



45 eV ^ a^ib eV)^^^ ^) 



. (0) 



■ (0) 



(0) 



1 - !(fy7) + (1 + !) (f^)(i + ^)^ - (fy.)(i + ^p 



(39) 



In Fig. [2] the 2;-dependent effective equation of state parameter Wejf is plotted for the 
case a = 2, P = 1 and various values of fh 



(0) 
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Figure 2. The plots show the effective dark energy equation of state parameter Weff 
as a function of redshift for different values of the neutrino mass scale rh^ — X]z=i '^i- 
The plots show the case a — 2, /3 = 1 and the left panel includes examples with 
ifil < 0.2 eV whereas in the right panel values ffii > 0.4 eV are presented. 



4-5. Comparison with Present Observations: Averaged Dark Energy Equation of State 

For a dark energy equation of state which varies with time, a model independent 
extraction of Weff{z) is difficult from the currently available data. Therefore, in most 
analyses of the present data, a constant w is assumed. One quantity which is used in 
the literature to compare with the best fit result under the assumption of a constant w 
is the weighted (or averaged) equation of state parameter Wavg defined in the following 
way [31] 

_ lo^effizWDEi^) dz 



W. 



avg 



j:n^j^liz)dz 



(40) 



.(0) 



In Fig. [3|, Wavg is plotted as a function of rh)^'' for the case a = 2,j3 = 1. The 
WMAP 5 year 2a limit on the constant dark energy equation of state parameter is 
—0.097 < 1 + w < 0.142 for zero curvature and —0.11 < 1 + w < 0.14 when a nonzero 
value of the curvature is allowed [31] • From the plot we can see that even ml as large 
as 1.5 eV is compatible with the present data. 



5. The MaVaN Hybrid Scenarios with High Scale Seesaw Mechanism 

One interesting result of the previous section is that with power law MaVaN potentials 
as in Eq. (jl]), the predictions for the observables do not depend on the value of A 
today. With the masses Mji of the right-handed neutrinos originating from terms in the 
potential of the form 

£mh = IxANN (41) 
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Figure 3. Averaged equation of state parameter Wavg as a function of neutrino mass 
scale fhi, — X]i=i '^j- The plot show tlie case where a = 2 and [3=1. 



this means that in principle the masses of the right-handed neutrinos today (~ AAq) 
could be close to the GUT scale and we might realise the MaVaN Hybrid Scenario with 
a "high-scale" seesaw mechanism. Compared to the usual MaVaN scenario where the 
right-handed neutrino masses are very small and where tiny neutrino Yukawa couplings 
have to be postulated, in the "high-scale" seesaw the neutrino Yukawa couplings can 
be of (9(1) and the smallness of the neutrino masses is explained by the large masses of 
the right-handed neutrinos. 

Recalling Eq. (120|) . the scales A^ and M in the considered scenarios are connected 
to the neutrino mass scale rhl by 



AnM 



^^^ I P - (0) (0) 



a 



l/a 



(42) 



To check the adiabaticity condition for the case of large Aq, we consider the square of 
the effective mass of the acceleron field today which is given by 



^lff = 



<fv,_ 



dA^ 



'^^ -a{a + (3){(3/af-^ 



p{o)i-^M§(4-a) 



Using Eq. (142|) in Eq. (l43l) we can eliminate the dependence on M and obtain 



m. 



eff 



A 



^(/3(a + /3)mWnW)i/2. 



(43) 



(44) 



Now, for adiabaticity to hold, rn^ff must be larger than the present value of the Hubble 
scale ~ 1.5 x lO^^^eV^. Assuming for example a = 2 and /? = 1 (i.e. a simple mass term 
potential for A and the standard seesaw relation), with Aq around the GUT scale ~ 10^^ 
GeV and rhi = 1.5 eV, rnejf is about 2 x 10"'^^ eV, which is considerably heavier than 
the Hubble scale. For the same choice of parameters and present vev of acceleron field, 
but with rhu = 0.06 eV; wig// ~ 4 x 10~^^ eV, still larger than the Hubble scale. This 
is in contrast to quintessence-type models, where the mass of the relevant scalar field 
has to be below the present Hubble scale. 
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We would however like to note at this point that although the scenario has several 
attractive features, we have not attempted the construction of a viable model of this 
type. One challenge is the radiative stability of the potential. Another challenge is the 
origin of the constant vacuum energy contribution. Both issues might be addressed in 
supersymmetric models (for example in simple setups of the form of Eq. ([6])), however 
even if hidden sector supersymmetry breaking is mediated to the SM fields by gauge 
interactions and to the dark sector only by gravity, this would typically induce a too 
large mass for A of at least about 10^'* eV. Leaving these model-building questions 
aside, we will now comment on how the "high-scale" seesaw might help with respect to 
stability issues for non-relativistic neutrinos. 

6. High Scale Seesaw and Stability Issues for Non-relativistic Neutrinos 

In this section we will comment on the stability issue regarding the neutrino density 
perturbations, which are typically considered as a serious problem for MaVaN scenarios 
with non-relativistic neutrinos. Let us first briefly review the problem: Soon after the 
introduction of the original MaVaN scenario [7], it has been pointed out [27] that the 
original scenario faces a catastrophic instability for non-relativistic neutrinos due to the 
extra force carried between neutrinos and the acceleron field A. With the mass m^ of 
the field A much larger than the Hubble scale, it has been argued that the adiabatic 
perturbations at scales between m^^ < A < H~^ are unstable. 

Let us now discuss how this problem might be overcome in the MaVaN Hybrid 
Scenarios with a "high-scale" seesaw mechanism. Rather than attempting a full 
numerical analysis of the density perturbations of the neutrinos in this paper (which 
we will leave for a future study), we will provide qualitative arguments and analytical 
estimates which will show how the MaVaN Hybrid Scenarios might help to suppress 
dangerous instabilities. We will start with a general consideration regarding the coupling 
between the neutrinos and the acceleron field, and then turn to analytical criteria for 
the occurrence of the instabilities. 

One main difference between the MaVaN Hybrid Scenarios with power law type 
potentials (where a "high-scale" seesaw mechanism is in principle possible) and the 
standard MaVaN scenario is that in the former the present value Aq of the acceleron 
field can be much larger than the eV scale, even close to the GUT or Planck scale. One 
can readily see that this might help with respect to the stability issue, when the heavy 
right-handed neutrino fields A^ are integrated out the theory. Then, the interaction term 
between the light neutrino(s) and the acceleron field becomes (for, e.g., /3 = 1) 

C-\nt = i^Li^L—^, (45) 

where A = Ao + 6 A with 6 A denoting the quantum fluctuations around the classical 
value of the field Aq. With 6 A <C Aq, the term in Eq. (l45l) takes the form 

-(0) 
171 

Ci^t^i^Li^L^SA, (46) 
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(C\\ 2 2 

where fhu = ^i^. In the usual MaVaN scenario the field value Aq is not much larger 



than ml ■ In contrast to this, as we have discussed in the previous section, in the MaVaN 
Hybrid Scenarios with power law potentials the present value of Aq can in principle be 
as large as the GUT scale (or Planck scale) and the coupling between neutrinos and A 
can be strongly suppressed (by a factor fhu / Aq). 

Let us now turn to the analytical criteria: First of all we note that when revisiting 
the stability issues of the MaVaN scenario, in Ref. [53] Bjaelde et al. have argued 
that a negative value of the sound speed does not always indicate the occurance of an 
instability. Due to the dragging force of cold dark matter, neutrino perturbations can 
remain stable even if the sound speed of the dark energy fiuid becomes negative. More 
specifically, the dragging force due to the dark matter and baryons can stabilise the 
neutrino perturbations if the condition ^5] , 

^CDM + ^b\ I G \ ( ScDM \ -. / . -x 



^u J \Geff 

is satisfied, where 

a,;^Gfl + ^ I. (4. 

V I + wi^"(A) + pJ'iA)) ' ' 

and where the coupling function f{A) is defined as 

fiA) = — ^ , (49) 

m^ dA 

with f{A) = —-J for the case /? = 1. Using the above definitions the stability condition 
of Eq. (HTl) can be rewritten as 

fiA) < C"'^':^ ?;,-'"' ] (^-^) . (50) 



2n,Ml, I \ 5, 



From this equation we can see that a small value of Vty (i.e. a small neutrino mass scale 
today) as well as a larger value of 5cdm compared to 5y can help to stabilise the neutrino 
perturbations. 

Setting ScdmI^u ~ 1 for simplicity and using Eq. fl50l) . we find that, e.g., for 
frbu ~ 1.5 eV, f{A) has to be smaller than 2Mp^ or for fhu ~ 0.06 eV, f{A) has to be 
less than lOMp^^. In order to avoid the instability the present value of A should therefore 
be larger than about Mp;/10 for ml ~ 0.06 eV and Mpi/2 for fhl ~ 1.5 eV. For the 
case ml ~ 0.06 eV this condition is well compatible with the adiabaticity condition in 



Eq. ( l22l) and the simple analytical consideration suggests that the instability problem 
might be cured. On the other hand, for fhl ~ 1.5 eV the adiabaticity and the stability 
conditions cannot safely be satisfied simultaneously and we conclude that whether the 



instability appears has to be checked numerically Jff 



ft We note that the calculation of sound speed (see e.g. 06]) does not depend on the present value of A 
and it is therefore still negative in our model, however, as we have mentioned before we are following 
the argument of j35j that despite a negative sound speed the instability can be avoided by the dragging 
force of dark matter. 
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Finally, we would like to note that, as we have discussed in section 4, for the most 
standard types of potentials (e.g. with a = 2, (3 = 1) the dynamical contribution to the 
total dark energy is rather small (typically a few percent), and the main contribution 
arises from a constant vacuum energy density. This means that even if the neutrinos 
would cluster, the observable effects on the smoothness of dark energy are much smaller 
than in the standard MaVaN scenario. Another related question is whether the scales on 
which neutrinos might cluster finally leads to any observable effect [37] . We conclude by 
remarking that although we have given some arguments which suggest that stability 
problems for non-relativistic neutrinos could be resolved in certain MaVaN Hybrid 
Scenarios with "high-scale" seesaw, this issue requires further investigations (which are 
left for future studies). 

7. Summary and Conclusions 

Motivated by the intriguing proximity of the energy density of dark energy and the 
neutrino mass scale we have studied the phenomenology of hybrid scenarios of neutrino 
dark energy, where in addition to a so-called Mass Varying Neutrino (MaVaN) sector, a 
cosmological constant (from a false vacuum) is driving the accelerated expansion of the 
universe today. Within the generalised framework we have focused on phenomenological 
issues such as on the connection to the neutrino mass scale and on its consequences for 
the dynamical nature of dark energy. 

We have therefore calculated the effective equation of state parameter Weff{z) in 
the MaVaN Hybrid Scenario where the effective potential for the dynamical real scalar 
field (the "acceleron" field A) has the following form 

We found that, for the case of a power law potential in the MaVaN sector, Weff{z) is 
determined by the neutrino mass scale and by the parameters a and (3 (c.f. Eq. ( I39l) ). 
Due to the interactions of the dark energy field with the neutrino sector, Weff{z) is 
predicted to become smaller than —1 for increasing 2; > (c.f. Fig. [21), which could be 
tested in future cosmological observations. For the considered scenarios, we have also 
calculated how the neutrino mass scale determines which fraction of the dark energy is 
dynamical, and which originates from the "cosmological constant like" vacuum energy. 
In particular, for the case of a mass term potential and a standard seesaw relation (i.e. 
a = 2 and /? = 1) we found that compatibility with the terrestrial neutrino mass bounds 
requires a large contribution of constant vacuum energy (c.f. Eq. (!38|) ). 

Another interesting question, which we have investigated in the MaVaN Hybrid 
Scenario with power law potentials is whether it is possible to realise neutrino dark 
energy with a "high-scale" seesaw mechanism, where the right-handed neutrino masses 
are close to the GUT scale. We found that the field value of the "acceleron" field as well 
as the masses of the right-handed neutrinos can indeed be large and in principle a hybrid 
scenario of neutrino dark energy might be realised with a "high-scale" seesaw. We have 



F(A),^y = M^-"A° + ^ -^ +K). (51) 
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also commented on how the Hybrid MaVaN Scenarios with "high-scale" seesaw might 
allow to suppress the formation of "neutrino nuggets" and resolve stability problems of 
dark energy models with non-relativistic neutrinos. 

In summary, we have found that the considered hybrid scenarios of neutrino dark 
energy have several attractive features, in particular the close connection to the neutrino 
mass scale and an effective dark energy equation of state parameter Weff{z) which 
depends only on the parameters a and (3 of the potential and on the neutrino mass 
scale. The prediction that Weff{z) < —1 for 2; > provides a "smoking gun" signal 
for such interacting dark energy scenarios, which could be observed in future surveys. 
Issues which are still open, and which are left for further studies, are whether a MaVaN 
Hybrid Scenario with "high-scale" seesaw can indeed solve the stability problems of 
conventional MaVaN scenarios with non-relativistic neutrinos and whether a consistent 
model can be constructed where this is realised. 

8. Acknowledgements 

We would like to thank Neal Weiner for useful discussions. This work was partially 
supported by The Cluster of Excellence for Fundamental Physics "Origin and Structure 
of the Universe" (Garching and Munich). SD is supported by US NSF CAREER 
grant PHY-0449818 and DOE OJI program under grant DE-FG02-06ER41417. SD 
acknowledges the hospitality of Max-Planck-Institut fiir Physik (Werner-Heisenberg- 
Institut) in Munich where this work was initiated. 

Appendix 

A. Some Remarks on the possible Origin of the MaVaN Hybrid Potential 

A.l. Fractional Power Potential from Non- Canonical Kinetic Energy Terms 

One possible origin of fractional powers in the potential in Eq. ([5]) is a non-canonically 
normalised real scalar field A' with the following Lagrangian 

A/n' 

^ = ]^ (5.^f - n^') , (A.l) 

where V{A') = Mf' A"^' ■ 

In the MaVaN Hybrid Scenario, the effective energy density of the field A' has the 
form 

yeff{A') = ^{d,Af + Mt-'A"^' + P-^[^^^ , (A.2) 

where the last term arises form the smooth background energy density. We have assumed 
m,{A') 



rNJ 



1 



A"^" 
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Now we would like to write it in terms of canonically normalised field A, which 
satisfies 

1 . . fA'-^"'!^ 



^S„A=^-j a^A'. (A.3) 

Integrating this equation we obtain the canonically normalised field A in terms A' , and 
the effective potential in terms of the normalised field A becomes, 

ln\ 2/3' 

n'a' „ / rS ) / A \ 77+2 

VeffiA) = g-'Mt-'^Mr-A^^ + ^ (t) ' ^^"'^ 



2 3 



where g = {n' + 2) "'+2 /2"'+2 which can be easily absorbed in the redefinition of M. 

As a consequence of the non-canonically normalised scalar field we have obtained 
fractional powers in the MaVaN Hybrid potential with a = ^^ and P = ^^r^- 

A. 2. Generalisation to Three Families: Possible Origin o/Vq 

One motivation for introducing a constant term Vq to the potential arises in 
super symmetric "hybrid-type" models of dark energy. Generalising the superpotential of 
Eq. ([6]) to three families, we may realise the situation that some of the "waterfall fields" 
Ni are still in the false vacuum while others are in the true vacuum where |A^j| ~ v^-. 
A simple superpotential with this characteristic feature is the following 

W = XA{x,N^ -v%) + VaiH^Ni . (A.5) 

In this case, the vacuum energy of the acceleron field does not account for the total dark 
energy, but only for part of it. The other part stems from the iVj which are in the false 
vacuum leading each to a contribution of energy density X'^vf^.. It is also interesting 
to note that the Ni which have non-zero vevs generate a mass term for A from the 
contribution |-F/vp to the scalar potential. 

A.3. Higher Powers of A from the Superpotential 



It is furthermore possible to generalise the setup of Eq. (JA.SP such that higher powers 



W = A— — (x.iVf -v')+ Vo^^HlKN, . (A.6) 



of A arise in the potential. Such higher powers can emerge from a superpotential of the 
form 

A'' 

The most relevant parts of the scalar potential (for a N^ in the true vacuum) would then 
be given by |-FaP ~ f^JAp""^ and l-Fjy.p ~ | ^jt-i Ni + yai-ff °z>q, p. This generahsation of 
the superpotential has two effects: Firstly, I-F4P leads to a higher power of A, explicitly 
|^|2n-2^ i.e. a = 2n — 2 in the effective potential of Eq. (jl]) and secondly, l-F/vJ^ results in 
masses of the right-handed neutrinos of the form m^, ~ ^\Jn-i which leads to masses 
of the light neutrinos fhy ~ ^^in v"^. This implies (3 = n m. Eq. (j4]). We note that 
another way to realise higher powers of (3 would be to generate neutrino masses via 
a multiple seesaw (see e.g. [30] for a model of neutrino dark energy using the double 
seesaw mechanism). 
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